Localized and Extended States in a Disordered Trap 
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We study Anderson localization in a disordered potential combined with an inhomogeneous trap. 
We show that the spectrum displays both localized and extended states, which coexist at intermedi- 
ate energies. In the region of coexistence, we find that the extended states result from confinement 
by the trap and are weakly affected by the disorder. Conversely, the localized states correspond to 
eigenstates of the disordered potential, which are only affected by the trap via an inhomogeneous 
energy shift. These results are relevant to disordered quantum gases and we propose a realistic 
scheme to observe the coexistence of localized and extended states in these systems. 
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Disorder underlies many fields in physics, such as elec- 
tronics, superfluid helium and optics [1-3]. It poses 
challenging questions, regarding quantum transport [4] 
and the interplay of disorder and interactions [5]. In 
this respect, ultracold gases offer exceptionally well con- 
trolled simulators for condensed-matter physics [6] and 
are particularly promising for disordered systems [7]. 
They recently allowed for the direct observation of 
one-dimensional (ID) Anderson localization of matter 
waves [8-11]. It should be noticed however that ultracold 
gases do not only mimic standard models of condensed- 
matter physics, but also raise new issues which require 
special analysis in its own right. For instance, they are 
most often confined in spatial traps, which has significant 
consequences. On the one hand, retrieving information 
about bulk properties requires specific algorithms [12]. 
On the other hand, trapping induces novel effects, such 
as the existence of Bose-Einstein condensates in low di- 
mensions [13], and suppression of quantum tunneling in 
periodic lattices [14]. 

Consider Anderson localization [15]. In homogeneous 
disorder, linear waves can localize owing to coherent mul- 
tiple scattering, with properties depending on the system 
dimension and the disorder strength [1]. A paradigm 
of Anderson localization is that localized and extended 
states generally do not coexist in energy. This relies on 
Mott's reductio ad ahsurdum [1]: Should there exist a lo- 
calized state and an extended state with infinitely close 
energies for a given configuration of disorder, an infinites- 
imal change of the configuration would hybridize them, 
forming two extended states. Hence, for a given energy, 
almost all states should be either localized or extended. 
Exceptions only appear for peculiar models of disorder 
with strong local symmetries [16]. Then, a question 
arises: Can inhomogeneous trapping modify this picture 
so that localized and extended states coexist in energy? 

In this Letter, we study localization in a disordered 
potential combined with an inhomogeneous trap. The 
central result of this work is the coexistence, at interme- 
diate energies, of two classes of eigenstates. The first class 



corresponds to states which spread over the full (energy- 
dependent) classically allowed region of the bare trap, 
and which we thus call "extended" . The second class cor- 
responds to states of width much smaller than the trap 
size, which are localized by the disorder, and which we 
thus call "localized" . We give numerical evidence of the 
coexistence of extended and localized states for different 
kinds of traps. We show that while the extended states 
are confined by the trap and weakly affected by the disor- 
der, the localized states correspond to eigenstates of the 
disordered potential, which are only affected by the trap 
via an inhomogeneous energy shift. Finally, we propose 
an experimentally-feasible scheme using energy-selective 
time-of-flight (TOF) techniques to observe this coexis- 
tence with ultracold Fermi gases. 

Let us consider a d-dimensional gas of noninteracting 
particles of mass m, confined into a spatial trap Ft(^) 
and subjected to a homogeneous disordered potential 
V(r) of zero average, amplitude and correlation length 
(Jr. Hereafter, we use "red-detuned" speckle potentials 
(Fr < 0), which are relevant to quantum gases [7, 17]. 
For the trap, we take VT{r) = {h^ / 2mo?)\r / ^ being a 
the trap length scale. For instance, a = oo and a = L/2 
for a homogeneous box of length L, while a = 2 and 
a = y/h/muj for a harmonic trap of angular frequency 
uj. We numerically compute the eigenstates \ipri) and 
eigenenergies En of the Hamiltonian 

H = -h^V'^/2m + V(r) + Vrr{r). (1) 

The eigenstates are characterized by their center of mass, 
rn = (V^nl^lV^n), and spatial extension (rms size), Ar^ = 

((^/^nl^^lV^n) — ''"n)^ ^- ^hc quantity Arn quantifies local- 
ization: the smaller, the more localized. 

Numerical results for the ID (d=l) case are reported 
in Fig. 1. In infinite, homogeneous disorder (a = oo, 
L = oo), all states |?/^n) are localized, uniformly dis- 
tributed in space, and, for most models of disorder, 
their extension Azn monotonically increases with the en- 
ergy [18]. As Figs. l(a),(f) show, using a finite flat box 
(L < oo) only induces a trivial finite-size effect: For low- 



Figure 1: Extension (a)-(e), center of mass (f)-(j) and density of states [DOS] (k)-(o) of the eigenstates versus energy in various 
kinds of ID disordered traps. The plots result from accumulation of numerical data over 5000 realizations of a red-detuned 
speckle potential with ttictr | Vr | //i^ = 0.256. The first column refers to a flat box of length L = 500crR. The curved line in (a) 
corresponds to an infinite system [28] and the horizontal line is Az*^ = L/2v^. The other columns refer to inhomogeneous traps 
with a = 12.5crR and various trap powers a. The solid lines correspond to the nondisordered case, i.e. Az^ in panels (b)-(e) and 
ztZci(E) in panels (g)-(j). The last row shows the full DOS p{E) (solid black line), as well as the DOS restricted to localized 
(p<, solid red line) and extended (y9>, dashed blue line) states [19]. The dot-dashed green lines are the nondisordered limits. 



enough energy we find Azn <C L and the states are 
not significantly affected by the finite size of the box. For 
larger energies however, boundary effects come into the 
picture. The states are centered close to the box center 
and their extension saturates to the value obtained for a 
plane wave, i.e. Az^ = L/2^/3. A central outcome of these 
results is that the curve giving Azn versus E displays a 
single branch. In particular, there is no energy window 
where localized and extended states coexist. This find- 
ing holds independently of the finite box size and is in 
agreement with Mott's argument [1]. 

For inhomogeneous traps {a < oo), we find a com- 
pletely different behavior. The curves giving Azn and 
Zn versus E now display two clearly separated branches 
[see Figs. 1(b)- (e) and (g)-(j)]. For low energy, the states 
are strongly localized and, for E > 0^ they are roughly 
uniformly distributed in a region bounded by the (energy- 
dependent) classical turning points, Zc\{E), defined as the 
solutions of Vt{zc\) = E. For higher energy, the exten- 
sion of the states corresponding to the upper branch in 
Figs. 1(b)- (e) grows and eventually saturates to that of 
the eigenstates of the nondisordered trap, Az^{E). The 
centers of mass of these states approach the trap center 
and form the horizontal branch in Figs. l(g)-(j). This 
branch corresponds to extended states. It is easily in- 
terpreted in terms of finite-size effects, similarly as for 
a finite fiat box. The lower branch in Figs. 1(b)- (e) is 
more surprising. It identifies strongly localized states 
of relatively large energy. It has no equivalent in the 
fiat box and cannot be interpreted as a finite-size effect. 
The corresponding states are located close to the classi- 
cal turning points Zc\{E) and generate the outer branches 
in Figs. l(g)-(j). As Fig. 1 shows, this holds for all in- 
homogeneous traps. When the trap power a increases, 
the branch of extended states gets denser at the expense 
of that of localized states, and completely vanishes for 
a = oo (homogeneous box). 

The coexistence of localized and extended states in the 
same energy window for disordered traps is confirmed on 
more quantitative grounds in the last row of Fig. 1. It 
shows the full density of states (solid black line), as well 
as the density of localized (p<, solid red line) and ex- 
tended (p>, dashed blue line) states [19]. The different 



nature of the localized and extended states is even more 
striking when one studies the wavefunctions. Let us fo- 
cus for instance on the harmonic trap (a = 2) and on a 
narrow slice of the spectrum around E ~ 4|]/r|, where 
p^/p 14% of the states are localized [20]. Figure 2(a) 
shows the spatial density |?/;n(^)P of all states found for 
a single realization of the disorder. We can clearly distin- 
guish localized (thick red lines) and extended (thin blue 
lines) states, which shows that they coexist in the same 
energy window for each realization of the disorder. The 
localized states are very narrow and present no node (e.g. 
states A and E) or a few nodes (e.g. states C and H). They 
may be identified as bound states of the local deep wells 
of the disordered potential, similarly as the lowest-energy 
states creating the Lifshits tail in bare disorder [18]. To 
confirm this, let us decompose the eigenstates of the 
disordered trap onto the basis of the eigenstates \xp) of 
the bare disordered potential [i.e. Hamiltonian (1) with 
Vt =0], associated to the eigenenergies Cp. For a local- 
ized state IV^n), we find KXplV^n)!^ ~ 1 for a single state 
\Xp) such that En, where = + (XpI'^t(^)IXp) is 
the eigenenergy of \xp) shifted by the trapping potential 
[see Fig. 2(b)]. Conversely, the same decomposition for 
an extended state shows a broad distribution of ampli- 
tude much smaller than unity. A localized state Itjjn) of 
the disordered trap thus corresponds to a strongly local- 
ized state \xp) in the bare disorder, which is affected by 
the trap by just the energy shift {Xp\^t{z)\Xp) - We gen- 
erally find that \ep\ <C {Xp\Vt{z)\Xp) - (V^n|V'T(^)|V^n), 
and, due to the reduced spatial extension of IV^n), we 
get En — Vrizn) [21]. This explains that the localized 
states are located close to the classical turning points, as 
observed in Figs. 1 and 2(a). 

Let us now decompose the states {t/jn) of the disor- 
dered trap onto the basis of the eigenstates li/j^) of the 
bare trap [i.e. Hamiltonian (1) with V = 0], associated 
to the eigenenergy E^. For a localized state, the distri- 
bution is broad. Conversely, for an extended state, the 
distribution is sharp and peaks at E^ En to a. value 
equal to a fraction of unity [see Fig. 2(c)]. An extended 
state may thus be seen as reminiscent of an eigenstate 
of the bare trap, which is weakly affected by the disor- 
der. Still, the main peak in Fig. 2(c) is smaller than 
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Figure 2: Eigenstates for a single realization of a ID disor- 
dered harmonic trap, (a) Non- normalized spatial densities, 
\ipn{z)\'^, vertically displaced to their eigenenergy En- Thick 
red lines correspond to localized states, and thin blue lines to 
extended states [19]. Note that extended and localized states 
may occupy almost- degenerate energy levels (e.g. H and I). 
The states C and D are projected: (b) over the eigenstates 
of the disordered potential, \xp)j and (c) over those of the 
harmonic trap, \ipp). The parameters are as in Fig. 1. 



unity. Only for significantly higher energy, the state {i/jn) 
results from weak perturbation of and |('0p|'0n)P 

displays a main peak of the order of unity as predicted 
by standard perturbation theory. 

Our results can now be easily interpreted. In bare dis- 
order, the typical size Az of the localized states increases 
faster than the classically allowed region Zc\ cx pro- 
vided by the trap. For low energy, Az <C z^ so that the 
states are strongly localized by the disorder and weakly 
affected by the trap. For higher energy however, the dis- 
order would localize the states on a scale exceeding Zc\. 
The states are then bounded by the trap and the effect 
of disorder becomes small. This forms the branch of ex- 
tended states in both the disordered box and traps. In 
addition, some strongly localized states with very low en- 
ergy in the bare disorder and located around point Zn are 
shifted by the trap to approximately the energy Vrizn)- 
This forms the branch of localized states only in disor- 
dered traps {a < oo) since, in the box, a state cannot 
be placed at intermediate energy due to the infinitely 
sharp edges. Quantitatively, since the localized states in 
the bare disorder are uniformly distributed in space, the 
density of localized states can be estimated to roughly 
scale as p< oc {I /a) x which is consistent with 

the disappearance of the branch of localized states when 
a grows and with its vanishing for a = oo (see Fig. 1). 
Still, it is striking that localized and extended states can 



coexist in the same energy window. The disordered po- 
tential combined with a smooth trap permits localized 
states to sit outside the classically allowed region occu- 
pied by extended states [see Fig. 2(a)]. Then, the Mott 
argument does not apply here because the spatial seg- 
regation can be strong enough to suppress hybridization 
for an infinitesimal change of the disorder configuration. 

Let us now discuss a possible scheme to observe the 
coexistence of localized and extended states in a disor- 
dered trap. Consider a gas of noninteracting ultracold 
fermions prepared in a given internal state, at temper- 
ature T and chemical potential fi. A class of energies 
\En — E\ < A [see Fig. 3(a)] deep in the Fermi sea 
(i.e. with (1 — E <^ k^T) can be selected by applying 
a spin-changing radio- frequency (rf) field of frequency 
V = E/h and duration r ^ h/ A (with h the Planck con- 
stant) [14, 22, 23]. The rf field transfers the atoms of 
corresponding energies to an internal state insensitive to 
the disordered trap. The transfered atoms expand freely, 
which provides their momentum distribution: 

VE,A{k) ^ E \Mk)\^ (2) 

\Er^-E\<A 

where t/Jnik) is the Fourier transform of V^n(^) [TOF tech- 
nique]. In the coexistence region, VE^/^ik) has two sig- 
nificantly different contributions: For localized states, 
|?/^n(/c)p is centered around k c:^ {) with tails of width 
Akn ~ Conversely, for extended states, \il)n{k)\^ 

is peaked k ^ \phnElh with long tails towards small 
momenta. We however found that averaging over real- 
izations of the disorder blurs the central peak associated 
to the localized states in VE,A{k). In turn, the quan- 
tity CEMk) = ^EMk) X P£;,a(0)/P£;,a(0)2 displays 
two distinct peaks for a rf pulse of realistic durations 
[see Fig. 3(b)]. The central one is more pronounced for 
narrower pulses. Selecting either the localized states or 
the extended states [19] confirms that the central peak 
corresponds to the localized states and the side peak to 
the extended states [see Inset of Fig. 3(b)]. 

Finally, we have performed similar calculations as 
above in a 2D harmonic trap. Figure 4(a) shows the 
centers of mass Tn of the eigenstates with En ^ 4|Fr|, 
the color scale giving Arn- Figure 4(b) shows a density 
plot of Arn versus \rn\ for the same data. Again, the 
eigenstates clearly separate into two classes: Some states 
are extended (large Arn) and centered nearby the trap 
center (small \rn\)- The other states are strongly local- 
ized (small Arn) and located nearby the line of classical 
turning points {\rn\ — ^ci(^) = \/ "IE j muj^\ Hence, the 
two classes of states can coexist at intermediate energies 
also in 2D disordered traps. 

In conclusion, we have shown that, in a disordered in- 
homogeneous trap, localized and extended states can co- 
exist in a given energy window. The localized states cor- 
respond to eigenstates of the disordered potential which 
are only affected by the trap via an inhomogeneous en- 
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Figure 3: Scheme to observe the coexistence of locahzed and 
extended states in disordered traps (sohd red hne). (a) Atoms 
occupying the eigenstates of energy EztA (shaded region) are 
transfered to a different internal state via rf couphng. The 
corresponding momentum distribution is then measured by 
TOF. (b) Correlation function CE,A{k) (black solid line) and 
momentum distribution VE,A{k) (dashed green line, arbitrary 
units), for A = 2h(jj. Inset: CE,A{k) of all states (solid black 
line), and separating localized (dashed red line) and extended 
(dotted blue line) states [19], for A = O.Ol/ia;. Here E = 41^1 
and the other parameters are as in Fig. 1. 
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Figure 4: Coexistence of localized and extended states in a 2D 
disordered harmonic trap for £^/|Vr| = 4 ±0.0003. The figure 
results from accumulation of data from 2 x 10^ realizations 
of the disorder, with mailVnl/fl^ = 0.8 and uj = 0.05\Vn\/h. 
(a) Centers of mass of the eigenstates and corresponding 
values of Arn/cTR in color scale. The solid black line is the line 
of classical turning points, rci{E) = ^2E/mu'^ ^ 63.2crR. (b) 
Extension Arn versus distance from the trap center \rn\- 

ergy shift. Conversely, the extended states spread over 
the classically allowed region of the trap and are weakly 
affected by the disorder. This effect is directly relevant to 
present day experiments with disordered quantum gases, 
which are most often created in harmonic traps [11, 24- 
27]. We have proposed a realistic scheme to observe it 
in these systems. In the future, it would be interesting 
to extend our results to higher dimensions and to other 
kinds of inhomogeneous disordered systems. 

We thank T. Giamarchi and B. van Tiggelen for dis- 
cussions and ANR (Contract No. ANR-08-blan-0016-01), 



Triangle de la Physique, LUMAT and IFRAF for sup- 
port. 



[1] N.F. Mott, Metal- Insulator Transitions (Taylor and 

Francis, London, 1990). 
[2] J.D. Reppy, J. Low Temp. Phys. 87, 205 (1992). 
[3] E. Akkermans and G. Montambaux, Mesoscopic Physics 
of Electrons and Photons (Cambridge University Press, 
New York, 2006). 
[4] A. Lagendijk et a/., Phys. Today 62, 24 (2009); A. As- 
pect, and M. Inguscio, Phys. Today 62, 30 (2009). 
[5] T. Giamarchi and H.J. Schulz, Phys. Rev. B 37, 325 
(1988); M.P.A. Fisher et a/., Phys. Rev. B 40, 546 (1989). 
[6] M. Lewenstein et al, Adv. Phys. 56, 243 (2007); 

I. Bloch et a/.. Rev. Mod. Phys. 80, 885 (2008). 
[7] L. Sanchez-Palencia and M. Lewenstein, Nature Phys. 6, 
87 (2010); L. Fallani, C. Fort, and M. Inguscio, Adv. At. 
Mol. Opt. Phys. 56, 119 (2008). 
[8] B. Damski et a/., Phys. Rev. Lett. 91, 080403 (2003). 
[9] L. Sanchez-Palencia et a/., Phys. Rev. Lett. 98, 210401 
(2007); New J. Phys. 10, 045019 (2008). 

[10] J. Billy et a/.. Nature (London) 453, 891 (2008). 

[11] G. Roati et al, Nature (London) 453, 895 (2008). 

[12] T.-L. Ho and Q. Zhou, Nature Phys. 6, 131 (2010). 

[13] D.S. Petrov et a/., J. Phys. IV (France) 116, 5 (2004). 

[14] L. Pezze et a/., Phys. Rev. Lett. 93, 120401 (2004); 
L. Viverit et a/., Phys. Rev. Lett. 93, 110401 (2004); 
H. Ott et a/., Phys. Rev. Lett. 93, 120407 (2004); 

[15] PW. Anderson, Phys. Rev. 109, 1492 (1958). 

[16] S. Kirkpatrick and T.P. Eggarter, Phys. Rev. B 6, 3598 
(1972); Y. Shapir et a/., Phys. Rev. Lett. 49, 486 (1982); 
O.N. Dorokhov, Solid State Comm. 51, 381 (1984). 

[17] We obtained the same qualitative behavior as reported in 
the Letter using other models of disorder, such as "blue- 
detuned" speckle potentials and Gaussian impurities. 

[18] I.M. Lifshits et a/.. Introduction to the Theory of Disor- 
dered Systems, (Wiley, New York, 1988). 

[19] We use the condition Az < IOctr to identify localized 
states, and Az>{3/4:)Az^{E) for extended states. 

[20] We obtained qualitatively-similar results for different en- 
ergies and different traps. 

[21] For instance, we find En ^ 3.971141, ^ 0.24|Vr| and 
{Xp\Vt\Xp) - 3.711141 for state C in Fig. 2. 

[22] M. Greiner et al, Nature (London) 426, 537 (2003). 

[23] W. Guerin et a/., Phys. Rev. Lett. 97, 200402 (2006). 

[24] D. Clement et a/., Phys. Rev. A 77, 033631 (2008). 

[25] Y.P Chen et a/., Phys. Rev. A 77, 033632 (2008). 

[26] L. Pezze et a/., Europhys. Lett. 88, 30009 (2009). 

[27] B. Deissler et al, Nature Phys. 6, 354 (2010). 

[28] Here, we use a box of length 3000crR much larger than 
Az for all states considered. 



